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In this work we study the large scale structure bias in models of anisotropic inflation. We use 
the Peak Background Splitting method in Excursion Set Theory to find the scale-dependent bias. 
We show that the amplitude of the bias is modified by a direction-dependent factor. In the specific 
anisotropic infiation model which we study, the scale-dependent bias vanishes at leading order when 
. . . the long wavelength mode in squeezed limit is aligned with the anisotropic direction in the sky. We 

CO ' also extend the scale-dependent bias formulation to the general situations with primordial anisotropy. 

, We find some selection rules indicating that some specific parts of a generic anisotropic bispectrum 

■ is picked up by the bias parameter. We argue that the anisotropic bias is mainly sourced by the 
angle between the anisotropic direction and the long wavelength mode in the squeezed limit. 

^ ■ I. INTRODUCTION 

(N : 

,_i . Inflation [ij has emerged as the loading paradigm for the theory of early universe and structure of formation. Basic 
' predictions of inflation indicate that the curvature perturbations are nearly scale- invariant, nearly adiabatic and 
, nearly Gaussian which are in very good agreements with cosmological observations such as WMAP [2| . The simplest 
models of inflation are based on a scalar field rolling slowly on a flat potential. On the other hand, the observations 
from Planck satellite are expected to improve the constraints on primordial non-Gaussianity(PNG) significantly. Any 
detection of primordial non-Gaussianity will have significant implications for inflationary model buildings, for a review 
O . see 0-0]. For example, many models of single field infiation predict a very small amount of local non-Gaussianity in 
the squeezed limit [6| Jnl ^ {^s — 1) in which Ug is the curvature perturbation power spectrum spectral index and 

■ /nl parametrizes the amplitude of local non-Gaussianity (NG). With Hs ~ 0.97 from WMAP Q, one expects Jnl < 1 
for conventional models of single field inflation. However, this expectation is violated if the system is not reached 
the attractor regime 0, 01 if allows for a non Bunch-Davies initial condition Furthermore, inflationary 



' models with large non-Gaussianity predicts different shapes for bispectrum. Therefore, any detection or otherwise of 

, large primordial non-Gaussianity with different shapes will go a long way to rule out many inflationary scenarios or 

\^ ' put constraints on model parameters. 

CO , The most suitable cosmological observation to constrain the primordial non-Gaussianity is CMB. This is because the 

' perturbations in the last scattering surface are in the linear regime and the fingerprints of non-Gaussianity are mainly 



CO . preserved [l^ll3|- However, recently the interests in LSS observations and their implications for non-Gaussianity are 
' boosted due the theoretical findings of scale-dependent bias In general, the distribution of baryonic matter in 
, the Universe, mainly clustered in galaxies and the clusters of galaxies, is the fundamental observable of LSS [isl [lij. 
. . ' The distribution of galaxies and clusters of galaxies can be studied by a) : Probability distribution function (PDF) of 
^ , structures (mass function of structures), and b): The correlation functions, power spectrum and even higher moments 
^> ' of distribution. 

^ , In the standard theories of structure formation, the Gaussianity assumption plays a crucial role in finding the distri- 
_C^j bution of structures via the primordial density contrast distribution 17 1. Accordingly, changing the initial condition 
from a Gaussian to non-Gaussian primordial density perturbations will change the PDF of structures. This change 
mainly shows itself in the tail of distribution function. Consequently, this effect manifests itself mainly in the statistics 
of the clusters of galaxies in high mass and high redshifts distributions [l8l - [2^ . 

The emergence of non-Gaussianity will introduce a non-zero bispectrum in LSS observations, however the gravita- 



tional instabilities are also the source of higher order nonzero moments j23j . The gravitational instability has its own 
signatures on the bispectrum of structures, but its shape changes in deep non-linear regime. So the bispectrum of 
galaxies and clusters of galaxies is not a very efficient observational tool to detect the primordial non-Gaussianity. 
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However, the primordial non-Gaussianity changes the background linear density fluctuations, which in turn it also 
influences the density peaks where structures are formed. The effect of background density fluctuations, affecting 
the halo abundance of dark matter, is known as bias parameter. The non-Gaussianity has a very unique feature in 
LSS by introducing a scale-dependent bias [l4, 24 1. This special fingerprint of primordial non-Gaussianity provides 
the opportunity to constrain primordial non-Gaussianity using the power spectrum of galaxies. Also there are other 
LSS observables such as the Integrated Sachs Wolfe cross correlation with the galaxy power spectrum [25|, the 3D 
bispectrum of Ly-alpha forest and the redshifted 21-cm signal from the post re-ionization epoch 26|, |27|, the statistics 
of voids [2^ etc. which can also be used to study the primordial non-Gaussianity. In a word, the LSS observations 
will become very important tools, complementary to CMB observations, to constrain the properties of primordial 
non-Gaussianity. 

There have been some indications of statistical anisotropics on CMB power spectrum. Although The statistical 
significance of the violation of statistical isotropy is not high [29l . [sot . but nonetheless the possibility of having 
statistically anisotropic seed perturbations are intriguing. The statistical anisotropy is usually parameterized via 
[sij Vc^ = 'PcqO- + g*{k.hY) in which n is the preferred direction in sky, Vc^ is the curvature perturbation power 
spectrum of the Fourier mode k with the direction along the unit vector k and V r^ r epresents the isotropic power 



spectrum. Constraints from CMB and large scale structure indicate that jgi*! < 0.4 3^ 



and the trispectrum in a model of anisotropic inflation [3J] have been calculated in [35 



3311 . Recently the bispectrum 
3g. It has been shown that 

large non-Gaussianities with non-trivial shapes are generated. Considering our motivation in using non-Gaussianity 
flngerprints in LSS as a probe of inflationary universe, we would like to study the effects of large scale-dependent and 
orientation-dependent bispectrum on bias in these models. 

The rest of the paper is organized as follows. In Section |TT] we present the anisotropic inflation model and the 
corresponding bispectrum which is used in subsequent analysis. In Section IlIII we review the basics of halo bias with 
non-Gaussian initial conditions. In Section HVl we present our results of halo bias in the models of anisotropic inflation 
with the anisotropic bispectrum obtained in Section |lll In Section |Vl we present a general mathematical formulation 
for the anisotropic bispectrum in terms of spherical harmonics which can be used to calculate the halo bias in models 
with generic anisotropic bispectrum. We leave some technical issues of halo bias analysis into appendix. 



II. ANISOTROPIC INFLATIONARY MODEL 



In this section we review the anisotropic inflation model and its anisotropic bispectrum which will be used in 
subsequent analysis. 

The best method to introduce statistical anisotropics is to incorporate U{1) gauge fields or vector fields in models 
of inflation. However, due to conformal invariance of the standard Maxwell theory in an expanding background, the 
background gauge field and its quantum fluctuations are diluted during inflation. Therefore, in order to produce an 
almost scale-invariant power spectrum of gauge field fluctuations, one has to consider a time-dependent gauge kinetic 
coupling. This prescription was originally used in 3^ 4^ in the context of primordial magnetic fleld. An interesting 



model of anisotropic inflation is introduced in |34l | in which it is shown that, with a suitably chosen gauge kinetic 
coupling, the inflationary system has an attractor solution in which the gauge field energy density reaches a small but 
obscrvationally detectable fraction of the total energy density. As a result, the anisotropy produced are at the order 
of slow-roll parameters. 

The Lagrangian of the system is 



(1) 



in which (f) is the inflaton field and F^^, = dfj^A^, — d,yA^ is the field strength associated with the U (1) gauge field A^. 
The background is in the form of Bianchi I universe with the metric 



= -df + a{tfdx^ + b{tf{dy^ + dz^). 



(2) 



Here H = a \s interpreted as the average Hubble expansion rate. Ha = a/a and Hb = b/b are the expansion rates 
along the spatial directions x and y and &/H = {Hb — Ha)/H is a measure of anisotropic expansion. We note that 
this metric enjoys only a two-dimensional rotational symmetry in y — z plane. 
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The details of the dynamics of the system are given in [34, 37|. For the simple chaotic potential V = m^(f)'^ /2, the 
conformal coupling (the time-dependent gauge kinetic coupling) is chosen to be 

with c a constant. With c ~ 1, one can check that the system admits the attractor solution in which the ratio of 
the gauge field energy density in the form of electric field energy density is a small and constant fraction of the total 
energy density. Defining the fraction of electric field energy density to the potential energy as 

n^lim^, (4) 



during the attractor regime one obtains 



c-l I 

R=^^H = ^eH, (5) 



in which en = —H/H^ is the slow-roll parameter and / : 
The power spectrum of the curvature perturbation, 



fc3 



(CkCk') = (2^)'^c(k)'5=^(k + k') , Pc^^nlk), (6) 
for the particular anisotropic infiation described above was calculated in 35-3^ 41-4^ which has the form 

Pc(fc) = Po(l + .9*(fc.ri)') , (7) 
where the anisotropy parameter is given by 

= -2UN{ki)N{k2) . (8) 

Here N{ki) represents the number of e-folding that the mode of interest k leaves the horizon. In our notation, the 
number of e-folding is counted backwards from the time of end of infiation by a{N) = a/exp(iV) so iV < 0. For 
example, N — Ncmb = —60 to solve the fiatness and the horizon problem. As a result, N{k) is calculated to be 

N{k) - Ncmb = In ( . (9) 



CMB 



in which kcMB represents the comoving mode which leaves the horizon at Ncmb = —60 e- folds before the end of 
infiation. To satisfy the observational constraints from CMB and large scale structure jg*] < 0.3, we require / < 10~^. 
The bispectrum of curvature perturbations, ^2, fca), is defined via 

(C(fcl)C(fe)C(fc3)) = (2^f <53(fcl+fc2+4)Sc(fcl,fc2J3)- (10) 

The Bispectrum for the model of anisotropic inflation was calculated using in-in formalism in 35 , 3^ and using SN 
formalism in 37[ with the result 

B^ikiMM) = 288/A^(fci)iV(fc2)A^(fc3) (c(fci,A?2)Po(fci)^o(fc2) + 2perm.) . (11) 

Here the anisotropic shape function C'{ki,k2) is defined as: 

C{ki, k2)^(l- (fci.n)' - (k2.nf + (ki.n) (fe.n) (fci.fc2)) , (12) 



where n is the specific anisotropic direction in the sky. Note that in Eq. pT|) . Po(^'i) represents the isotropic power 
spectrum so all anisotropics are encoded in shape function C(fci,fc2) (and the appropriate permutations) with the 
amplitude 28&IN{ki)N{k2)N{k3). 
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It is instructive to look into the bispectrum in the squeezed hmit in which one mode is much longer than the other 
two. say fcs <C fci ~ k2, so from the condition ki ~ wc also conclude that ki ~ — fc2- In this limit, one obtains 



fc2,fc3) ^ 24Poik,)Poik3)\g,{ki)\N{k3) x 1 - ik,.hf - {k^.nf + {ki.n){k3.n){ki.k3) (fcg < k 



(13) 



in which to obtain the above result, Eq. ([8]) have been used to express the parameter IN{ki)N{k2) in terms of 
The non-Gaussianty parameter /jvl is defined in the squeezed limit fc3 <C fci ~ ^2 via 



/Ari(A:i, fc2, kz) = lim 



5 Bc_{ki,k2, ks) 

kT-^^O 12 P<;{ki)P^{k3) 



(14) 



In general, /nl is an orientation-dependent and scale-dependent quantity. As an order of magnitude estimation, and 
neglecting the logarithmic scale-dependence in N{ki), we can define an orientation-dependent effective f^/ via 



NL 



240/iV, 



CMB 



(fci)C(A;i,fe), 



60, we can easily get f^/ 



(15) 

60 with ~ 0.1, compatible 



where N{k) = Ncmb + ln(A:/fccA/s)- Setting Ncmb 
with observational constraints. 

A very interesting observation from Eq. (|13p is that when the long wavelength mode k^ is in the direction of 
anisotropy, i.e. k^ \\ n, then the term inside the big-bracket in Ea. (|13p vanishes. Consequently, in this configuration, 
we do not expect to see the non-Gaussian effects in LSS. We discuss this feature in more details in Sec lIIII 

For the subsequent analysis we adopt the coordinate system such that the anisotropic direction fi coincides with 
the z direction in the spherical coordinates so the other momentum vectors are described by fci = {Oijipi) and 
^2 — {02, 11^2), where 9 and ip are the polar and azimuthal angles in spherical coordinates, respectively. For the future 
reference, we also need the angles between two arbitrary unit vectors qi = {Oq., ipq.) defined via cos 7 = 171.(72, which is 



cos 7 = sin(0gj) sin{0q2) cos(-0gj — ipq^) + cos Oq-^ cos 



^q2 ■ 



(16) 



The power spectrum and the bi-spectrum presented in Eqs. ([7]) and (jlip are for the particular model of anisotropic 
inflation as studied in 35 -3^. For a generic anisotropic model the most general power spectrum can be written as 

M 



P(fc)=Po(fc) 



LM 



9LM{k)YLM{k) 



(17) 



where Pq is the isotropic power spectrum, Y^Mik) (with L > 2) are spherical harmonics and gLM{k) quantify the 
departure from statistical isotropy as a function of wavenumber k. As each Fourier mode k is related to — fc, in the 
case of real gLAiik), the multipole moment L must be even, and in the limit of fc — > we will recover the isotropic 
power spectrum Po{k). However, in the general case, (real and imaginary gLAd), we have: 



9lm (-1)^.9L-A/ . 



(18) 



This condition is imposed by the fact that the matter power spectrum is a real quantity. 

Comparing Eq. (|17p with Eq. ([7]) for our particular anisotropic inflation model we hvae 520 oc while the rest of 
9lm are zero. In Section |V] we extend the general deflnition of Eq. ()17p for the power spectrum to the bispectrum 
and look into its implications in halo bias analysis. 



III. BIAS 



In this section we review the concept of the bias, a parameter that shows the dependence of dark matter halo 
abundance to the background dark matter density perturbations. The reader who is familiar with these analysis can 
directly jump to next Section in which we present our results of halo bias for anisotropic inflation model. It is worth 
to mention that in this work we are not interested in galaxy bias, which is the weighted integral of the halo bias, 
corresponding to the mechanism of the halo occupation distribution (HOD). 
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We follow the work by Scoccimarro et al. [48|, to find an expression for the bias. The halo bias relates the halo 
abundance to the dark matter overdensity, in Excursion Set Theory (EST) , (also known as Extended Press Schcchter 
(EPS)) Hi], defined by: 

b{k,z) = ^, (19) 

where is the halo overdensity and 5m is the matter density perturbation. The EST framework is very useful as it 
provides a non-perturbative model (that includes e.g. exclusion effects at small scales )[5^. At large scales EST is 
known to reproduce the Gaussian initial condition while in small scales it determines the local bias parameter with 
linear and nonlinear terms [sil [H^ which are in reasonably good agreements with numerical simulations (HsI - fssj . 

It is worth to mention that the halo bias is a function of redshift and scale. This scale-dependance is introduced 
by applying the initial Non-Gaussian condition. Bias parameter in the context of peak-background splitting (PBS) 
sil . 53] is described by the fact that the background large scale overdensity changes the critical threshold of spherical 



collapse [58| such that the collapse criteria will be : 

Ss>5c-Si, (20) 

where Ss is the matter density contrast of the structure, (the subscript "s" stands for short wavelength); Si is the 
background (long- wavelength) density contrast and 6c — 1.68 is the critical density contrast in spherical collapse 
formalism. 

Now in order to find the bias parameter we have to compare the dark matter halo abundance, in cases with and 
without the presence of long wavelength (background) overdensity. For this task we use the Excursion Set theory 
approach. In the Appendix [XJ we review the concept of EPS in more details and we will derive the bias parameter 
using the PBS in EPS context. For a review of PBS sec Appendix 1X1 

According to Appendix [A] the large scale halo biasing can be treated in peak-background splitting ^] . The idea 
of splitting of the density contrast to short and long wavelength can be translated to a similar splitting of Bardeen 
potential due to Poisson equation, which depends on cosmological parameters. The matter density in PBS can be 
written as: 

p^pil + 6s + 6i) . (21) 

The number density of formed structures with mass m can be expressed as a function of small scale statistics, (i.e. 
smalls scale power spectrum Ps{k)) and the background long wave-length perturbation Si (i.e. n = n[Si{x), Ps{ks)\ to]) 

The primordial potential in Fourier space can be translated to the late time potential as: 

^k,z) = ^,„,xT{k)Diz){l + z), (22) 

Where <f>i„i is the initial Bardeen potential sourced by the inflation field quantum fluctuations, T(k) is the transfer 
function and D{z) is the growth function normalized to scale factor at early times. Now due to Poisson equation one 
can relate the initial non-Gaussian potential to Si via Poisson equation in sub-horizon scale and linear-regime: 

Si=M{k,z)^, (23) 

where 

, , 2k^Tik)D(z) , , 

Here, fio and i?o are the matter fraction energy density and Hubble parameter at present time, respectively. 

In order to calculate the halo-bias term, we should calculate the effect of large scale perturbation. Si. on the 
probability distribution function(PDF) of structures. This yields the following relation between the Lagrangian halo 
number density and PDFs 



^ ^ ^ n(,5,, g^, Suaf)dS, ^^^^ 

^™ /oJ no {Ss , cr?^ , S^dSs 

where, H((5s, cr^, (5c; <5;, f;^) is the conditional PDF, which means that the variance Om at large scales converges to the 
value Gra = c/, in contrast with unconditional PDF, Ho((5s, cr,^„, (5c), in which the variance vanishes (om ^ 0) at large 
scales. 
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It will be relevant to define a quantity that shows the probability of first up-crossing in the time interval between 
(T^ and (7^ + dcr^j, in EST language as: 



d 



So 



ToiS,,ai) = -j— Ilo{S,,ai,dc)dS,, (26) 

O'^ni J-oo 

From the above formalism, we can see the effect of primordial NG on LSS. Assuming that the Bardeen Potential 
has a local-type non-Gaussianity we have 

$ = + /JVL0' (27) 

where </> is the Gaussian field. Now we can use the splitting idea on (f> by applying (/) = + (j)s, where (f>i is long wave- 
length mode of potential and (fi^ is short wavelength corresponding to the scale of structure. In Appendix \X\ we discuss 
how the above non-linear form can be generalized to a model with arbitrary shape of non-Gindependentaussianity, in 
which case, the non-linear term generalizes to a kernel. 

Now in the presence of primordial non-Gaussianity, the modes are not independent and the conditional PDF is 
modified by the non-Gaussian long-wavelength mode. In this case the PDF will be a function of (pi through the 
variance and also higher order cumulants (cp = {S^)c) as: 

Il{Ss,al,S,; Si,0)^ n[Ss,a^ {((>), Cpi4>),S,; 5,(0), 0], (28) 

That is, unlike the case of Gaussian initial condition that only the variance depends on the long- wavelength, in the 
case of non-Gaussianity the higher cumulants are also affected. Under the assumption that all these effects are small, 
and also using the Markovian random walk in EST formalism we can taylor expand the conditional PDF, 11, around 
unconditional one, Hq. For details, see appendix \X[ 

In this work we study the linear bias effect of anisotropic primordial non-Gaussianity and include only the first 
derivative contribution in taylor expansion, Ea. (|A15[) . As we show in appendix 1X1 up to first order in Jnl-, only the 
first two terms in taylor expansion of PDF appear. These terms are derivatives of PDF with respect to the long 
wavelength mode 6i and the variance cr;. The higher order terms corresponding to derivatives with respect to Cp 
{p > 3) contribute to 0{flf]^) and 0{gML) bias. 

As we mentioned, in this work we only study the first order effect in /tvl ■ The p — 1 contribution, the first term in 
taylor expansion, Ea. (jA16p . is the usual scale-independent linear bias from Gaussian perturbations. Keeping in mind 
h = Sfi/Si, for the first order linear bias (6il) we have 



p=l ■■ biL 

which can be written as: 



d^mJidU/dSi)o 

d,n J IIo 



_9_ j^^dn(^^ 
d6i d In m 



(29) 



biL = ^hi{n{Si)). (30) 
dSi 

In the presence of primordial non-Gaussianity, there are new contributions from higher order cumulants p > 2. As 
a result, the next to leading order term, gives 

p = 2. b2L = TTTTTt — Trf ' ^"^^^ 

M(k)drn J Ho 

which in general is a scale-dependent correction to the leading order, scale-independent bias (|A18[) . The key quantity 
here is J21 which is the derivative of second cumulant <t,^j, {p = 2), with respect to the long wavelength mode 0; which 
is obtained as [isj : 

hlik, m) = J B^^^iq, k-q, -k)d\, (32) 

where B^~^^ is the cross bispectrum of small-scale smoothed density 5 and 0. As a result, I21 is the quantity which 
we arc looking for in the case of non-Gaussian initial condition which introduces scale-dependent bias of the order of 
OUnl). 
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So far only the Lagrangian bias appeared because peaks are those of the initial density field (linearly extrapo- 
lated). Making the standard assumptions that halos move coherently with the underlying dark matter, and using the 
techniques outlined in one can obtain the final Eulerian bias as: 



bs = I + biL + b2L; (33) 
Now Assuming the Markovian behavior of density fields and universality of mass function we will have: 

2 

biL = ydina^^ ln(cr™ J") = 6il(g) + biL{NG), (34) 

A very important point is that in above equation we omit the subscript of T, which means that the non-Gaussianity 
will change the PDF of structures, so the first linear term will have a contribution from primordial non-Gaussianity. 
Consequently, we have: 

blL{G) = ^5ln<T^ M<^m^o) (35) 



and 

OlL(WG) = — "in^T^^ ln(cr„j7ijvG) = , (36) 



where = TZngJ'o, and TZng comes from the deviation of PDF from the Gaussian case [21|,|64|. In other words, the 
effects of non-Gaussianity appeared both in the PDF and the power spectrum via scale-dependent bias parameter. 
Since in this work we arc interested in the scale-dependence features of bias, the contribution of TZng is not much of 
interest. For the Gaussian case we use the Sheth and Tormen [6^ Gaussian mass function. For the non-Gaussian PDF 
effect we use the results of [l^ h^ which they used the idea of expanding the non-Gaussian PDF in the Press-Schechter 
framework [l7j in which TZng is defined as: 

TZNGim, Inl) = 1 + lx{x' - 3)ss{x) - ^(.t - l/x)^^, (37) 
o D aln(x) 

where x = Sc/tm and Sc = 1.68 is the critical density and S3 is the reduced skewness defined as: 

^3(^) = - (38) 

The skewness is related to the matter bispectrum as: 

(4)= / |^^M/(i?'7i)W^(i?<Z2)T^(i?gi2)M(gi,z)Af(g2,z)M(gi2,z)i?o(gi, 92,912), (39) 

where qi2 = Qi + Q2, W is the window function in Fourier space and R is the smoothing scale. For more detail see 
Appendix [A] 

On the other hand, the scale-dependent bias can be rewritten by 

/2i(fc,m) 1 /2i(fc,m) 

So the total Eulerian bias up to first order in /nl can be split into the scale-independent bsi and scale-dependent b^d 
terms as: 

bt(E) = bsi + bsd (41) 

where hsi and bsd are: 

bsi^ba + biUNG), (42) 

with 6g = 1 + 6iL(G) and 

bsd = b2L + 0{fj,L)+0{gNL) + - (43) 

In the next section we find the scale-independent and scale-dependent bias and discuss the effect of the primordial 
anisotropy in the bias parameter. 
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IV. ANISOTROPIC BIAS 



111 this section we study the LSS bias in the model of anisotropic inflation. For a related work with a phenomeno- 
logical modeling of bispectrum and its implication for bias see [66[ . 

In this work we assume that the mechanism of spherical collapse for structure formation is still applicable in our 
model with primordial anisotropic perturbations. This is motivated from the fact that the collapse mechanism is a 
local process so up to first order it is not affected by large scale anisotropics. Secondly, we assume that the transfer 
function at leading order is not affected by the anisotropic terms in energy-momentum tensor of cosmic fluid. It would 
be interesting to perform the analysis with the above assumptions are relaxed but this is beyond the scope of this 
work. 

For this task we investigate the change in PDF and corresponding cumulants, and the effect of these changes on 
bias. First of all the variance is defined through the linear matter power spectrum: 

^P^{k)W'{kR), (44) 

where Pl and W{kR) are linear matter power spectrum and window function in Fourier space which arc discussed 
the details in Appendix [Xj 

In our model the variance is modified due to the fact that we use anisotropic power spectrum. However, we show 
below that it is not direction-dependent. Without loss of generality, we can assume that the anisotropy is pointed 
along the z-direction in spherical coordinates {h — z). Then, starting from Eq. ([7]) for the primordial anisotropic 
power spectrum, we have 

cr^(m, z) = J di^dcosBdkk^Poik, z)[l + g^{k) cos^ 9]WikR) ~ (1 + y)crf„ (45) 

where 6 and ip are spherical coordinate angles and (7*(fc) = —2AIN{ki)N{k2)- Note that (j\ is the variance obtained 
from the full anisotropic power spectrum whereas (t^„ is the variance corresponding to the isotropic part. Since the 
scale-dependence of 5* is logarithmic through N{ki), as given in Eq. ©, we can ignore its scale-dependence as a first 
approximation, so we have the last approximate equality in Ea. (|45p . Note that, as we mentioned before, the variance 
does not have any direction-dependence, which is somewhat an obvious observation, since one should integrate over 
the full 3D Fourier space to obtain the variance, eliminating any direction present in power spectrum. However, it 
is interesting to note that the correction to the variance due to anisotropy enhances or suppresses the leading order 
term, depending on the sign of anisotropy parameter . From the above variance, we can obtain the leading order 
scale- independent bias 

As for the next step, we obtain the bias from the anisotropic bispectrum which is both scale-dependent and direction- 
dependent. This would be the main effect of anisotropy in the bias parameter. In order to obtain 62L, we need I21 
where the information from primordial bispectrum is encoded. We use the squeezed limit (fca ^ ^2 ~ fci) bispectrum 
predicted by the model, Eq. p^ . Because of symmetry in x — y plane we rotate the long wavelength mode such that 
'0^ = 0. Inserting the bispectrum to /21, Ea. ([32)) results in 



/2i(k,z,m) = 24jd^gd{cos0g)dqq'M„^iq,z)Mrn{\k-q\,z)P^{q)N{k)\g,{q)\ (46) 

X [1 — cos^ 9q — cos^ 9k + COS 9q COS 9k (sin 9q sin 9k cos -0, -I- cos 9q cos 9k)\ 

where k and q correspond to long and short wavelength (k^ and k2 in Eq. ()13p ). and 9, ip are polar and azimuthal 
angles in spherical coordinates, defined by the angle between the anisotropy direction h and wavenumbers q and k 
respectively: 

cos7q = q.h = sm9n sm9q cos{tpn — ipq) + cos 6n cos 9q, (47) 

cosjk = k.h = sm9nS'm9f,cos{ipn — ip^) + cos9nCOs9f,, (48) 
Now we can integrate the angular dependence jq which yields 

/2i(A:,z,m,0fc) = 647riV(fe)(l-cos2 0fc) f dqq^ MmAl)MmA\k - q\)P4'iQ)\9*i<l)\ ■ (49) 
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Anisotropic NG, Q=90 deg 

Anisotropic NG, Q=7 deg 




k (h/Mpc) 

FIG. 1: The fraction of Scale dependent anisotropic bias Ea. (|55|l . to the Gaussian bias, b[,'^'' /ba, is plotted versus the wavenumber 
for different angles between the long wavelength mode and the anisotropy direction. The black double-dashed curve, the green 
long-dashed curve and the blue dotted curve, respectively, are for 6k ~ 33° , 6^ = 90° and 9k ~ 7° . The solid red line indicates 



the local non-Gaussianity with /jvl = 100. For the anisotropic bias, we choose I 
(fT5l) . The redshift is 2 = 0. 



1.9 X 10"" such that /^^ = 100 from Eq, 



Anisotropic NG, Q=33 deg ■ 
Local NC.f„i=100 ■ 



k (IVMpc) 



FIG. 2: This figure is a zoom-in of Fig. [T]for long wavelength 
modes where we compare b^^V^G of our model with 6k = 33° 
and the local non-Gaussianity with /nl ~ 100. This shows 
small deviation between the two models due to a mild scale- 
dependence of parameter. 



Local NG, f,j^im ■ 
Anisotropic NG, 0-.?J deg - 



k (h/Mpc) 



FIG. 3: This figure is a zoom-in of Fig. [T]for short wavelength 
modes where again we compare fe^^V^G of our model for dk ~ 
33° with the local non-Gaussianity with /nl ~ 100. 



It is interesting to note that the orientation-dependence appears in the form of sin^ 9k. As a result, the bias vanishes 
when sin^^fe = 0, i.e. when the long wavelength mode is aligned with the anisotropic direction. This result originates 
from the fact that in this specific direction the bispectrum vanishes in squeezed limit as one can check from Eq. (jl3p . 
Furthermore, since in the model under consideration has a mild scale-dependence via logarithmic correction in 
N{k), we observe that there is an extra, but mild k-dependent, factor in I21 in comparison with the standard local 
non-Gaussian shape (48l l67j . Besides that, I21 linearly depends on / which is the free parameter of the anisotropic 
inflationary model. Now, by using the variance and I21 parameter obtained above we can find the bias parameter in 
the anisotropic model by 



dm /Ho 



d6i d In m 



^IL(G) + biL{NG): 



(50) 
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Q=33 cleg, 1=10 
Q=33 deg. 1=5*10'' 
e=ii deg, 1=10'' 



T 1 1 1 1 — I I I 

6 




—I I I I I ' ' ' ' 



_l I I l_ 



_l I I I I I— l_L 



k (h/Mpc) 

FIG. 4: &i;^VbG is plotted versus wavenumber with 9^ = 33° for three different values of /. The solid red curve, the green 
dashed curve and blue dotted curve, respectively, are for I = 10~^ / = 5 X 10"'' and I = 10"^ In all case we set the redshift 
z=0. 



— 1 1 1 1 — I — I I I I 



— 1 1 1 — I — I I I I — 

Q=33deg, z=0 
Q=33deg, z=l 
Q=33deg, z=2 



T 1 1 1 1 — I I I 



—I I I I I ' ' ' ' 



_l I I l_ 



_l I I I I I— l_L 



k (h/Mpc) 

FIG. 5: feld'/^G is plotted versus wavenumber with 7 = 10 ^ and 9k = 33° for three differnt redshifts. The solid red curve, the 
green dashed curve and the dotted blue curve, respectively, are for z — 0, z — 1 and z — 2. 



where, for the first order linear bias in the case of anisotropy, we have 

, 2 2t:, n dlnTZ(^A)NGimjNL) 

0{A)1L{NG) = g-^^'^A ^^'y'^A''^{A)NG) = 



dSi 



(51) 



in which the subscript (A) has been added to point out that the parameters are obtained in the presence of anisotropy. 
Here TI{^a)ng is the anisotropic non-Gaussian correction to the PDF of primordial fluctuations defined by: 



'^{A)NGi'm, fNh) = 1 + \xa[x\ - 3)s(a)3(xa) - \{XA - 1/xa)^^4~^' 

6 ^ ' D aln(a;^j 

where xa = 5c/cta and 5c = 1-68 is the critical density and 5(^^)3 is the reduced skewness defined as: 

(4) 



S{A)3iR) = T?2 

\"{A)bJ 



/2 



(52) 



(53) 
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Note that, in the above formula, we have to use the full variance a a including the correction due to anisotropy. The 
anisotropic skcwness is also defined as: 

{^Ur)= j -^-^WiRqi)W{Rq2)W{Rqi2)M{qi,z)Miq2,z)M{qi2,z)BA{qi,q2,qi2), (54) 

where Ba is the anisotropic bispectrum. 

The next order term in bias will give the scale-dependent and direction-dependent effect: 

AA) , /, n\ l2l{k,Z,m,e) 1 /2l(fc,TO,6>fc) 

b,, =M)2L(fc,^,™,^.)=2^^(^^^^),^(^^^)'5AL + ^g(^a,„..( ) (55) 

So the total Eulerian anisotropic bias is defined as: 

bt{E) = 1 + biL(G) + b(A)lLiNG) + b(A)2L, (56) 

The first three terms above give the scale-independent bias, whereas the last term is scale-dependent and direction- 
dependent bias due to the anisotropic bispectrum. Since the model we consider has a bispectrum shape very close to 
the standard local non-Gaussian shape and since &iL(ArG) is small in comparison with bsd in local shape, |67| . we can 
ignore 6ii(jvG) term in our analysis. 

In Fig. [T]we plot the relative magnitude of the scale-dependent bias Eg. ([55]) to Gaussian bias, 6|^^^^/6g, versus 
the wavenumber. As expected we have approximately fc~^ scale-dependence similar to the local non-Gaussian shape. 
However, this fc-dependence is slightly different because of mild dependence of to wavenumber. A completely new 
feature is the direction-dependence of bias originated form primordial anisotropy, proportional to sin^ ■ In the case 
of ~ 7r/2 we have the maximum scale-dependent bias, whereas at angles = 0, tt, as the anisotropic bispectrum 
vanishes, the scale-dependent bias also vanishes accordingly. 

In order to compare our results with the conventional local non-Gaussian models, we set Ok^ = 33° to have 
~ 100 and then compare the bias in this specific direction with the bias in local shape with amplitude Jnl ~ 100. 
In Fig. Oand Fig. [31 for both long and short wavelengths, we compare the fraction of scale-dependent bias, b(scL)/bGi 
from local non-Gaussianity with our anisotropic model. Since oc N{k) = Nqmb +^'n-{k/kcAiB), for long wavelength 
modes the bias of our model is slightly higher than local non-Gaussian case while it is slightly lower for short wavelength 
modes. 

The free parameter of our model is /. In Fig. 2] we plot b^^^^^/bQ versus wave number for different values of /. The 

/-dependence of bias is linear according to Eq. In Fig. [S]we investigate the redshift-dependence of b'^^^-^/bc- 

At higher redshifts we have higher bias. This is the standard expectation, since at earlier times the non-linearity of 
local, small scale perturbations were weaker and thus they are more sensitive to large scale perturbations, resulting 
in larger bias. 

V. GENERAL FORMULATION 

In this section we calculate the bias, assuming a general model independent primordial anisotropy. For this task 
we write the power spectrum and bispectrum in their most general anisotropic forms. Then we find the possible 
configurations that give the scale-dependent bias. 



As mentioned in section |TI1 the general anisotropic power spectrum can be written as 47 



P{k) = Po{k) 



l + J29LMik)YLMik) 



LAI 



(57) 



which is a generalization of our anisotropic power spectrum defined in Eq.(I7]), where all the fc-dependence is buried 
in the coefficient gLM- The anisotropic variance reads as: 



LM 



(58) 



crl(m,z)= J d^d{cose)k^dkPo{k,z)W{kR) 
The spherical harmonics are orthonormal via 

J YimYi*^,dn ^ Sll'Smm' ■ (59) 
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As a result, noting Yqq = l/\/47r, we have 

(60) 

so the variance, Eq. ([55|) . can be simphfied to 

al{m,z) = 4:TT J k^dkPo{k,z)WikR) l+goo{k)/\^ . (61) 
Ignoring the scale-dependence of goo, we have 

(7\{m,z) ~ al{m,z) l+5oo/%/47r , (62) 

with 

cr^(m, z)^4:n J k^dkPo{k, z)W{kR) . (63) 
Similarly, we can generalize the bispcctrum in squeezed limit, where /ca <C fci w fc2, as follows 

B^{kiMM) ~ 2 CLlmVL{cOSe3)Yi.rni01,^Pl) (Po (^1 )/'o (^3 )) (64) 



where Vl is the Legendre function, Yi^miSi, 4'i) ^'^^ the spherical harmonics, by which the whole direction-dependence 
of bispcctrum is encoded and the cum are the fc-dependent coefficients. This decomposition is possible because we 
set the z coordinate along the anisotropy direction and rotate the x — y plane in Fourier (k-space) such that the 
azimuthal angle of k^ is set to zero. We can do this rotation because of the symmetry in the x — y plane. This kind of 
decomposition, i.e. decomposing into a fc-dependent sector(cLim) and the angular-dependent part is especially useful 
to find the corresponding I21 term in the case of generic anisotropic bispcctrum and the corresponding scale-dependent 
bias. 

Since the in the real space the bispcctrum is a real quantity, and noting that the bispcctrum above is defined in 
Fourier space, we find the condition 

ck„ = (-l)'^+W™. (65) 

Now, using Eq. (j64[) . we can write I21 as 

/21=2 ^ /^^l^(c0^^l)9''^'?^W'(9)AW(|fc-9l)cLim('7,fc)7'L(^3)i1,m(^l,V'l)i"o(g), (66) 

which, again using the orthonormality of spherical harmonics, simplifies to 

/2i(k,m,z) = 4\/^^7'L(cos6'fc) X / q^dqM,nA<lWmA\k - q\)cL,i=o,m=a{q,k)Po{q) . (67) 

L 

This an interesting result showing that the anisotropic bispectrum induces an anisotropic, scale-dependent bias, 
through the Legendre function of the angle between the anisotropic direction n and the long wavelength mode. It is 
also interesting to note that there is a Selection Rule for scale-dependent bias at 0{fNL), where the harmonic numbers 
of the short wavelength must vanish (^ = 0, m = 0). All the other parts of anisotropic bispcctrum with l,m ^ 0, do 
not contribute to the leading order scale-dependent bias. 

It is worth to mention that the reality assumption of the bispcctrum implies 

cloo = (-l)''CLoo, (68) 

where in the case of real cloo we have to choose even Ls. Now, if wc define a direction-independent parameter I^lpi 
as 

lft^2iik,m,z) = iV^ X j q^dqMmA<lWmA\k-q\)cLm{q,k)PQ{q), (69) 
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then I21 will be re-expressed as 

l2i{k,m,z,ek) = Y,T^Licos0k)I^l)2iik,m,z) (70) 

L 

where 0^. is the angle between h and k. Now the scale-dependent bias becomes: 

J2L'PL{cos0k)lfl)2iik,m,z) 1 EL'^i(cOS7)^(L)2l(^.™>^) 
&2L(fc, Z, m, Ok) = r-TTT^ — r SMl + Y7n — \^^nal ( r? ^ ) C^l) 



where in this case a\ is defined as in Eq. (|6T1) 

Our specific model of anisotropic inflation studied in previous Section corresponds toL = 2,/ = 0, m = with the 
corresponding coefficients. 



VI. CONCLUSION AND DISCUSSIONS 



The large scale structure observations, like the statistics of rare objects and the scale-dependence of bias parameter, 
can be used as complementary cosmological observations to CMB observations to constrain the inflationary models. In 
this work we obtained the scale and direction-dependent bias of dark matter halos in anisotropic inflationary models. 

The anisotropic model used in this work has a bispectrum shape very close to the local non-Gaussianity shape with 
an extra mild k-dependence and also a direction-dependence. We showed that the bias parameter is mainly influenced 
by the angle {Ok-^) between the anisotropy direction and the long wavelength mode in the squeezed limit by a 
factor sin^ ■ An interesting observation is that the scale-dependent bias vanishes in linear order along the direction 
Oks = 0, TT. This can be explained by the fact that the bispectrum vanishes in squeezed limit when the long wavelength 
mode is aligned to the anisotropy direction. This means that, at the level of bispectrum, the model reduces to the 
Gaussian model along this direction. 

As it is clear from the PBS formalism, the bias parameter is mainly affected by the long wavelength mode. We can 
see this explicitly in our work where the angle between the anisotropy direction and the long mode appears in the 
bias formula. 

The free parameter of our model is /. Recent CMB observations imply that / < 10^^. Roughly speaking our bias 
parameter analysis also imply the same order of magnitude for /. Considering the fact that the observations are now 
in good agreement with a Gaussian bias we expect that in quasi-linear regime (where we have the strict constraint on 

bias) the ratio of /be is at the order of one so from our results we also find / < 10~^. 

In this work we also formulated the general, model independent anisotropic scale dependent bias in linear order. 
For this purpose, by assuming rotational symmetry in a; — y plane of Fourier space, we modeled the bispectrum as 
a function of spherical harmonics, Yim/ Logendre functions, Vl, of the angles between the short/long wavelength 
modes and the anisotropic direction. Interestingly, we find a selection rule for scale-dependent anisotropic bias, which 
shows that the bias parameter only responses to the Yoo part of the anisotropic bispectrum. We also show that bsd is 
independent of short wavelength mode and its direction. This is an obvious check for the formalism, since the bias 
should be only a function of the long wavelength mode. 

In this work we assumed that the mechanism of spherical collapse and the conventional form of the transfer 
function are still applicable. This is motivated from the fact that the collapse of structures is a local mechanism 
which is less affected by the large scale anisotropics, although in principle it is an interesting question to see how one 
can generalize the spherical collapse mechanism in the presence of primordial anisotropic fluctuations. On the other 
hand, the modification of the transfer function due to Boltzmann and perturbed Einstein equations in the presence 
of a generic anisotropic cosmic fiuid is an intersting question which is beyond the scope of this work. 

Finally, it is an interesting open question to use the formalism we presented in this paper to compare more rigorously 
the predictions of models of anisotropic inflation to data from LSS observations to constrain the free parameters of 
these models. 
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Appendix A: Halos bias in the Peak Background Splitting in the context of Excursion Set Theory 



In this appendix, first we review the Excursion Set theory (EST), then we will apply the Peak Background splitting 
(PBS) in EST context and finally we argue that how we can find the halo bias. In the context of Excursion Set theory 
(EST) halo formation can be described as the random walk of matter density contrast as the smoothing radius goes 
from very large radius, (corresponding to infinitesimal variance ct^, small S) to the scales crossing the linear threshold 
for collapse 6c at some finite smoothing radius. This radius is related to the scale in which the halos form. Within 
the EST formulation the number density of collapsed objects (dark matter halos of mass m ) per unit mass is given 
by: 



dn 
dm 



Urr 



T\-Q{5s-,CF'^n^5c)d5s 



(Al) 



where cr^j is the variance of the small scale density field smoothed with filter (window function) at spatial scale R 
and p is the background energy density, relating the halo mass to smoothing radius by m = AirpR^ /i. Furthermore, 
no((5s, cr^„, (5c) is the unconditional probability distribution function (PDF) of small scale perturbations reaching 5s 
(short wavelength density contrast) at variance . By the subscript as well as the unconditional for PDF, we mean 
that the initial condition (first step in random walk) is (5^ = where cr^j ~ 0, and it satisfies the absorbing barrier 
condition TIq{5c, cr^, 5c) = 0. 

The smoothing procedure will be done by top-hat window function (Hsj 



W{x) 



3(sina; — a;cosx) 



The lagrangian radius of collapsed objects is related to the mass via: 



R = 



1.162 X W^K^MqVlI 



1/3 



Mpc 



(A2) 



(A3) 



where m is the mass of the structure and = ri^(l + z)^. 

In the Extended Press-Schechter(EPS) formalism the variance am is a monotonic function of mass scale so we use 
it as the clock. So it will be relevant to define a quantity that shows the probability of first up-crossing in the time 
ai and al, + dal,, as: 



Toi5c,a,^) 



2 \ — 



d 

m — oo 



^oiSs,al„5c)d5s, 



Now the number density of collapsed objects will be: 



/ dn 
\dm 



^\^\xTo{5c,al), 
TO dm 



(A4) 



(A5) 



It is worth to mention that the number density of dark matter halos obeys the normalization condition: 

/ (4^)mdTO = p (A6) 
J dm 

An important parameter in mass function of structures is the density contrast power spectrum, defined by [HsI 

Pi(fc) = Ak''''T^{k)D^{z), (A7) 

where Hs is the spectral index, A is the linear matter primordial power spectrum amplitude in k ~ 
T(k) the transfer function, and D(z) is the growth factor normalized to scale factor at early times. 

The evolution of density contrast is imprinted in growth function Di^z) and the transfer function T{k), showing 
the scale-dependence of gravitational potential during cosmic evolution. In this work we use the growth function of 
standard ACDM cosmology and the transfer function of Bardecn, Bond, Kaiser and Szalay(BBKS) [H^] respectively 
as below: 



D{z) = - 



5 1 



21 + z 



f^f,r-J7A + (l + ^)(l + ^) 



(A8) 
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and 

T{k qil^hHlpc-^) « ^"^^^'^g^^^^'^^ X [1 + 3.89g + {I6.2qf + {5A7qf + (6.71g)*] "'^^ . (A9) 

A crucial point to indicate is that the variance has a Unear dependence to growth function ct^(M, z) = a^{M)D^{z), 
where ct^(M) is the present value of variance. This redshift-dependance is important in the sense that the statistics 
of structure and even the bias parameter will be redshift-dependent. In upcoming section we will show the redshift- 
dependence of bias parameter. In the matter dominated era D{z) scales like 1/(1 + z) which is a decreasing function 
with respect to redshift. Now that we set the matter density variance and large scale statistics of matter distribution, 
by setting the initial condition of perturbations we can find the mass function of structures in the Universe. For the 
Gaussian initial condition the probability distribution function (PDF) in Fourier-space top-hat filter is: 

noiS,,al,S,)=rG{Ss,<yl)-VGi2S,~Ss,al), ( AlO) 

where Vg is the Gaussian PDF. 

The assumption of the universality of mass-function (i.e. Press-Schechter) yields 

no(<5„a^,5c) = n — ,— ), (All) 

which leads us to write c^^o just as a function of threshold quantity u — 5c/<7rn'' 

alMSc,<rl)^'-^, (A12) 



where /(;/) is the usual Gaussian factor in Press-Schechter theory jl7|. Now in order to find the halo bias term in 
Gaussian and non-Gaussian inflationary models, first wc discuss the PBS. We will show how this method is applicable 
for non-Gaussian fields. 

For inflationary models with primordial non-Gaussianity, the gravitational potential is usually written as [HqI - ItH 

^x) = 0(x) + InlMx)]' - {4>^)) (A13) 

where is a Gaussian random field. The above relation is the simplest extension of Bardccn potential to include non- 
linearity, called local non-Gaussianity. We c an g eneralize the local type non-Gaussianity to arbitrary shape, simply 
by replacing the non-linear term by a kernel j48| 

$(x) = 0(x) + fNLK[^{x),(j){x)] . (A14) 

In order to obtain the bias parameter, we can use the Peak-Background Splitting idea on the Guassian potential 
(j) = (ps + 4>h where <j)i is the long- wavelength mode of the potential and is the short-wavelength mode. Since 
the large scale perturbations change the background for small scale modes, the PDF of small scale perturbations is 
modified to conditional PDF, n((5s, cr^j, 5c\5i, af). In conditional PDF, the initial condition is (Jm cr; in large scale 
limit (first step in random walk), in contrast with the unconditional case in which the variance vanishes on large 
scales. 

On the other hand, the existence of any type of non-Gaussianity modifies the PDF such that it is no longer 
independent of higher cumulants, while in Gaussian case the PDF is completely determined by zeroth {5s) and first 
(ct^) order cumulants. 

By the splitting procedure explained above, and noting that the PDF is now a function of all cumulants, we can 
expand the Lagrangian halo ovcrdcnsity as an expansion over large-scale modes (isj , 

r r'- ^ATT 2 x y ^' + 2 / r ' ^ r'^ ^xu (A 2 a ^ '^'(fci)</';(fc2) + ... (A15) 

where (...)o means that the corresponding quantity is evaluated at (pi = 0. Note that because the background is now 
modified by large scale perturbation (pi, the cumulants are now functions of (pi. The first derivative in Eq. (jA15[) . to 
all orders in primordial non-Gaussianity, is 
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where the cumulants are defined by 

cW^^z, c^^^^al, c(^)EEa^(^0, ct^^{5l),, c^^^ ^ (S^^i))^ (A17) 

in which the subscript m for cumulants show that they are evaluated without the background 0/ (independent oi 4>i). 
It is worth to indicate that to first order in fpfLi only the first two terms in the taylor expansion above (p = 1,2) 
contributes to the bias, while p = 3 contributes to 0{ffjj^) and 0{gNL). 

Now we can obtain the bias parameter, using the above formulation. The p = 1 contribution is the usual scale- 
independent bias presented in initially Gaussian case. Keeping in mind b = 5h/ 5u we have 



Om J Ho 

which can be simplifies as 



d dn{Si) 
dSi d In m 



(A18) 



&iL = 4-ln(n(<50) (A19) 
oSi 

The p — 2 contribution is the scale-dependent correction to the leading order bias coming from primordial non- 
Gaussianity 

P ~ Z . b2L - . w, rTf (,A2Uj 

M(k)drn J Ho 

The quantity I21 includes the information about the primordial non-Gaussianity, and it is the derivative of second 



cumulant tr^, (p=2), with respect to long wavelength mode 0; which is obtained as 48 1: 

hlik, m) = J B^^^iq, k-q, -k)d\, (A21) 

where B~^~^^ is cross bispectrum of small-scale smoothed density 5 and (p. 

So far the Lagrangian bias appeared, because peaks are those of the initial density field (linearly extrapolated). 
By the standard assumptions that halos move coherently with the underlying dark matter, one can obtain the final 
Eulerian bias as: 

bE = l + hiL + b2L , (A22) 

using the techniques outlined in 60l-63l|. 

Note that, due to the existence of primordial non-Gaussianity, the leading order scale- independent bias also modifies. 
Assuming the Markovian behavior of density fields and universality of mass function we have 



2 

blL = ^^IntT^^ ln((Tm-7^) = 6lL(G) + biL{NG), {A23) 



where in the Ea. (|A23p we omit the subscript of which means that the Non-Gaussianity will change the PDF of 
structures, resulting in a modification of 6il. As a result we have: 

2 9 

biL{G) = ^5in<Ti ln(ff™^o) (A24) 

and 

blL{NG} = ^C)in^2^ ln(CT47^7VG) = , (A25) 

as described in Section 

TlNGimjNL) = 1 + lx{x' - 3)ss{x) - lix ~~ l/x)^^, (A26) 
6 6 aln(a;) 
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where x = Sc/ctau = 1-68 is the critical density and S3 is the reduced skcwness defined as: 



The skewness is related to the matter bispcctrum as: 



(5|)3/2 



(4) = / ^^W{Rq,)W{Rq2)WiRqi2)Miqi,z)M{q2,z)M{qi2,z)Bo{qi,q2,qi2), 



(A27) 



(A28) 



where qi2 = qi + q2 and W{kR) is the window function in Fourier space, smoothing perturbations up to scale R. 

In obtaining the PDF of non-Gaussianity model in this approximation, we have assumed that all the deviation is 
imprinted in the skewness which may not be entirely true. In order to improve the results, numerical simulations are 
done [22, T^, Z3|- Consequently, a scaling parameter k defined by RNcix) — >■ RNoit^x) are introduced where, in the 
work of [74j. from simulation of (75| . it is obtained to be k = 0.91. (For a similar correction from simulation you can 
see [li) 

On the other hand, the scale-dependent bias can be rewritten by 



_ l2i{k,m) 



/2i(fc,m) 

■Cln alS 5 j 



(A29) 



So the total Eulerian bias up to first order in fj^L can be splitted to scale-independent hsi and scale-dependent hsd 
terms as: 



'Jsi -I" Osd 



(A30) 



where 6,,- and are: 



with 5(3 = 1 + 6]^ 



L(G) 



and 



52L 



0(/ 



nl) 



0{g 



NL) 



(A31) 



(A32) 
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